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. . ^ Abstract. Millisecond and binary pulsars are the most stable astronomical standards of fre- 

0^ ■ quency. They can be applied to solving a number of problems in astronomy and time-keeping 

' metrology including the search for a stochastic gravitational wave background in the early uni- 

, verse, testing general relativity, and establishing a new time-scale. The full exploration of pulsar 

• properties requires that proper unbiased estimates of spin and orbital parameters of the pulsar 
O i' be obtained. These estimates depend essentially on the random noise components present in 
D , pulsar timing residuals. The instrumental white noise has predictable statistical properties and 

makes no harm for interpretation of timing observations, while the astrophysical/geophysical 
low-frequency noise corrupts them, thus, reducing the quality of tests of general relativity and 
decreasing the stability of the pulsar time scale. 
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c/^ : 

• Timing observations of single and, especially, binary millisecond pulsars are widely 
recognized as extremely important for progression in a number of branches of modern 

O ' astronomy and time-keeping metrology. In particular, the implication of pulsar timing for 
^ \ testing general relativity in the strong-field regime 1 and setting up the upper limit on 
^ ■ the energy density of the stochastic gravitational wave background in the early universe 
I I \ are among the most stimulating on-going research activities. Another important aspect 

of pulsar timing is metrology of time on long time intervals providing a time-scale based 
on the high-stable rotation of millisecond pulsars around its own axis ^i4j and/or center 
of mass of a binary system [5l [6] . 
^ The accuracy of pulsar timing observations is now approaching 10 nanoseconds. Such 
> high precision requires construction of adequate data processing software taking proper 
\ account of the relevant physical effects that can contribute to the variations of timing 
^\ ' residuals [7]. Usually, the procedure for estimating pulsar parameters is based on the 
' premise that white noise is a dominant source of instability in times of pulse's arrival. 
However, long-term monitoring of pulsars certainly reveals the presence of a non-white 
component of the noise having the external origin [8]. This noise is called red as it has 
a spectrum that diverges at zero frequency. The lower the timing activity of the pulsar, 
the further toward low frequencies one must look in order to detect the red noise in the 
timing residuals. The red noise also makes the residuals correlated on long time intervals. 
_ Developing a computer code accounting for effects of the red noise on timing residuals 
and estimates of parameters is worthwhile. It requires further improving the model of the 
red noise [9] and the methods of estimation of the stability of the pulsar rotational/orbital 
phase [im [TT] . The present state-of-the art statistical analysis of pulsar timing data in 
the presence of a red noise has not yet reached the required level of completeness and a 
more elaborate technique has to be invented. 

Analysis of the stability of pulsar's rotational and orbital phase is more informative in 
time rather than in frequency domain. This is because any noise contains both stationary 
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and non-stationary parts, while spectral analysis of noise in frequency domain is adequate 
only if the noise is stationary. A non-stationary part of noise affects observed values of 
pulsar's parameters, increases their variances, gives rise to larger correlations and makes 
some of the parameters biased [TO] . 

Another problem relates to the ergodicity of the overall timing process. There is only 
one, available to us, realization of the pulsar timing series. It means that one can not 
access the statistical ensemble to process the data. Instead, we have to apply time in- 
tegration to evaluate the parameters. If the ergodicity is violated the true value of the 
parameters will be biased, thus, affecting our ability to maintain the pulsar time scale 
and to test general relativity. 

Third problem is related to the mathematical structure of the estimators of pulsar's 
parameters. The estimators are integrals in the frequency domain with a kernel being the 
spectrum of the noise. Spectrum of the red noise is divergent, and so are the estimators. 
We do not know what is the best mathematical technique for making the red-noise 
estimators unbiased and convergent. Definite progress in this area was achieved within 
the theory of distributions [Hj . 

Continuous monitoring of pulsar rotational/orbital phase is not sufficient for establish- 
ing a pulsar time scale - an ensemble of pulsars, called a pulsar timing array (PTA), is 
necessary |12] . Currently, two PTAs are continuously monitored by Australian [13j and 
Russian [M] pulsar groups. A PTA can be also used as a gravitational-wave telescope 
that is sensitive to radiation at nanohertz frequencies, and for improving the planetary 
ephemeris |T5]. 
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